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Abstract: 

We present numerical results for total cross sections and various distributions for 
e + e~ — > WW — > 4/(+7) at a future 500 GeV linear collider, obtained from the Monte 
Carlo generator RacoonWW. This generator is the first one that includes 0(a) elec- 
troweak radiative corrections in the double-pole approximation completely. Owing to 
their large size the corrections are of great phenomenological importance. 
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One of the most important class of processes to be investigated at future e + e~ linear 
colliders are the four-fermion production processes, which involve in particular W-pair 
production, e + e~ — > WW — > 4f . In order to match the experimental accuracy for this 
process of 1% or better, theoretical predictions for the corresponding cross sections of 
e + e~ — > 4/ with a precision at the level of some 0.1% are needed. This requires to include 
the complete set of lowest-order diagrams for e + e~ — > 4/ and the 0(a) corrections to the 
W-pair production channels e + e~ — > WW — > 4f. In addition, the leading higher-order 
corrections should be taken into account. 

In all regions of phase space where W-pair production dominates the cross section for 
e + e~ — > 4f, an expansion of the matrix element about the poles of the resonant W prop- 
agators provides a reasonable approximation for the radiative corrections. Neglecting 
corrections to the non-doubly-resonant contributions leads to uncertainties of the order 
a/ ii x r w /M w x l°g(. . .) ~ 0.1% with respect to the leading lowest-order contributions, 
where the logarithm indicates possible logarithmic enhancements. This so-called double- 
pole approximation (DPA) provides a gauge-invariant answer and allows us to use the 
existing results for on-shell W-pair production [ III and W-boson decay [§,£§], as far as 
the virtual corrections are concerned. 

The DPA for 0(a) corrections to pair production of unstable particles has already 
been proposed in Ref. [|5|. Recently different variants of the DPA have been used in 
the literature [ §, [7|, ||, ^ |i~0fl . A Monte Carlo generator including the corrections to the 
W-pair production subprocess and the leading-logarithmic corrections to the W-boson 
decays has been constructed in Refs. [ 0, U, but non-factorizable corrections and W-spin 
correlations have been neglected there. With this generator also results for the centre-of- 
mass (CM) energy ^/s = 500 GeV have been produced. A different, but also approximate, 
implementation [ Q of the DPA has been compared at energies ^/s > 500 GeV with a high- 
energy approximation [ |ll|] for e + e~ — > WW — > 4/. The first complete calculation of the 
0(a) corrections for off-shell W-pair production, including a numerical study of leptonic 
final states for LEP2 energies, was presented in Ref. [ || using a semi-analytical approach, 
which is, however, only applicable to ideal theoretical situations. 

Recently the first Monte-Carlo generator that incorporates the complete 0(a) cor- 
rections for off-shell W-pair production in DPA has been constructed. This generator, 
called Racoon WW [ |TD , [T^|, includes the complete lowest-order matrix elements for 



e + e~ — > 4/ for any four-fermion final state. For the virtual corrections a DPA is used 
without any additional approximations. In particular, the exact four-fermion phase space 
is used throughout. The virtual corrections consist of factorizable and non-factorizable 
contributions. The former are the ones that are associated to either W-pair production 
or W-boson decay; the results of Refs. [ [Lj ||] are used in this part. The latter comprise 
all corrections in which the subprocesses production and decays do not proceed indepen- 
dently. Up to some simple supplements, the virtual non-factorizable corrections can be 
read off from the literature [ fEIJ Q ; we made use of the results of Ref. [ Q and included 
the exact off-shell Coulomb singularity [ |T5]. The real bremsstrahlung corrections are 



based on the full matrix-element calculation for e + e~ — > 4/7 described in Ref. [ [L6| . More 
precisely, the minimal gauge-invariant subset including all doubly-resonant contributions 
of the processes e + e~ — > WW — > 4/7, i.e. the photon radiation from the CC11 subset, 
are included. By using the exact matrix elements for the real radiation, we avoid prob- 
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lems in denning a DPA for semi-soft photons (E^ ~ Tw) and, moreover, can include the 
leading logarithmic corrections to non-doubly-resonant diagrams (background diagrams) 
exactly. The real corrections and the virtual corrections are matched in such a way that 
all infrared singularities cancel exactly. The initial-state collinear singularities are regu- 
larized by retaining a finite electron mass and factorized into lowest-order matrix element 
and splitting functions. The collinear singularities connected to final-state radiation are 
treated inclusively, i.e. photons within collinear cones around the final-state fermions are 
integrated over, so that no logarithmic final-state fermion mass dependence survives. All 
contributions have been implemented in two programs, one of which uses the subtraction 
method described in Ref . [ |17| , the other one uses phase-space slicing. All parts of the 
calculations have been performed in two independent ways. A detailed description of the 
calculation and the Monte Carlo generator RacoonWW will be published elsewhere [ 
12 1 . Results for the LEP2-energy region have already been given in Ref. [ pL0|| . Here we 



concentrate on results for a/s = 500 GeV, which is relevant for a future linear collider. 

For the numerical evaluation we used the same input as in Ref. [ |Kj, i.e. the fixed- 
width scheme and the following parameters: 



— 


1.16637 x lO^GeV -2 , 


a = 


1/137.0359895, 


M w = 


80.35 GeV, 


r w = 


2.08699 . . . GeV, 


M z = 


91.1867 GeV, 


r z = 


2.49471 GeV, 


m t = 


174.17 GeV, 


m h = 


150 GeV, 


m e = 


510.99907 keV. 







(1) 

The reliability of the fixed-width scheme at high energies has been demonstrated in 
Ref. [ [TjJ. The weak mixing angle is fixed by c w = Mw/Mz, = 1 — c^. The parameters 
in ([]]) are over-complete but self-consistent. Instead of a we use G^ to parameterize the 
lowest-order matrix element, i.e. we use the effective coupling 
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in the lowest-order matrix element. This parameterization has the advantage that all 
higher-order contributions associated with the running of the electromagnetic coupling 
and the leading universal two-loop m t -dependent corrections to the dominant contribu- 
tions are correctly taken into account. In the relative 0{ot) corrections, on the other hand, 
we use a, since in the real corrections the scale of the real photon is zero. The W-boson 
width given above is calculated including the electroweak and QCD one-loop corrections 
with a s = 0.119. We do not include QCD corrections to the process e + e~ — > WW — > 4/, 
and initial-state radiation is only taken into account in 0(a). 

In Table [l] we present numbers for total cross sections without any cuts based on 
20 million events. In particular, we give the CC03 cross sections, i.e. the cross sections 
resulting from the signal diagrams only (defined in the 't Hooft-Feynman gauge). We also 
give numbers resulting from the complete set of diagrams for those final states where this is 
possible without cuts, i.e. the CC11 class of processes. In the considered cases, the effects 
of the background diagrams are below 0.6%. Note that we treat the external fermions as 
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final state 


CC03 Born 


full Born 


CC03 corrected 


full corrected 




87.12(4) 


86.66(3) 


90.27(4) 


89.81(4) 




87.11(4) 


86.65(3) 


90.19(5) 


89.73(5) 




261.4(1) 


260.0(1) 


270.6(1) 


269.3(1) 




261.3(1) 


260.0(1) 


270.3(2) 


269.0(2) 


udsc 


784.1(3) 


780.3(3) 


811.3(4) 


807.5(4) 




784.0(3) 


780.4(3) 


810.9(4) 


807.2(4) 


total 


7057(2) 




7305(2) 






7056(2) 




7299(3) 





Table 1: Total cross sections in fb for e + e — > WW — > 4/ without cuts for various final 
states at y/s = 500 GeV 

massless. Therefore, the cross sections for processes not in CC11 class become singular if 
no cuts are imposed for final-state electrons collinear to the beams and for virtual photons 
splitting into // pairs with small invariant masses. The shown corrections are typically 
3.5% at \/s = 500 GeV. The numbers in parentheses are estimates of the Monte Carlo 
integration errors. The numbers for the total W-pair (CC03) production cross section in 
the last two rows of the tables directly result from the other rows by multiplying these 
with the number of equivalent channels and adding them up. The numbers in the upper 
lines are from the subtraction-method branch of RacoonWW, the number in the lower 
lines from the phase-space slicing branch. The agreement between these numbers serves 
as a stringent consistency check of our generator. 

Note that the total CC03 cross section of Table |l] is different from the full cross sec- 
tion for four-fermion production, as non-doubly-resonant (background) diagrams signifi- 
cantly contribute to the total cross section for high energies. The most important class 
of background diagrams are the ones for single-W production, which is dominant if elec- 



trons or positrons are scattered in the very forward directions (see e.g. Ref. [ [H| and 
references therein). The DPA approach is, of course, restricted to situations where the 
doubly-resonant diagrams dominate the cross section. Therefore, before confronting the 
(9(a)-corrected prediction of RacoonWW with empirical data one has either to specify 
cuts that suppress background diagrams or to extract the W-pair signal from the full 
four-fermion signature. 

Next, we study observables that involve phase-space and photon-recombination cuts. 
Because of our treatment of mass singularities (see Ref. [ |l!| for details) it is necessary 
to combine photons that are collinear to the incoming or outgoing fermions appropriately 
with these fermions in order to obtain well-defined finite distributions. To this end we 
introduce the following recombination and cut procedure which proceeds in three steps: 

1. All photons within a cone of 5 degrees around the beams are treated as invisible, 
i.e. their momenta are disregarded when calculating angles, energies, and invariant 
masses. 



3 



2. Next, the invariant masses of the photon with each of the charged final-state fermions 
are calculated. If the smallest one is smaller than M rec , the photon is combined with 
the corresponding fermion, i.e. the momenta of the photon and the fermion are added 
and associated with the momentum of the fermion, and the photon is discarded. 

3. Finally, all events are discarded in which one of the final-state fermions is within a 
cone of 10 degrees around the beams. No other cuts are applied. 

We consider the cases of a tight recombination cut M rec = 5 GeV, of a loose recombination 
cut M roc = 25 GeV, and the case of full photon recombination, M rec = 500 GeV. In the 
following observables, the momenta of the W bosons are always defined by the sum of the 
momenta of the two corresponding decay fermions after the eventual recombination with 
the photon. 

We now discuss various distributions for the z/ M /i + du final state. These results have 
been obtained from 50 million events. In the following figures we always show on the 
left-hand side the absolute distributions in lowest order and including the corrections 
for the recombination cut M rec = 5 GeV, and on the right-hand side the corresponding 
relative corrections for the three recombination cuts M rec = 5 GeV, M rec = 25 GeV, and 
M rec = 500 GeV. 

In Figure [l] we present the distribution of events in the angle between the W + and 
the incoming e + . The corrections lead to a drastic increase of the cross section in the 
backward direction. This effect is due to hard-photon emission from the initial state, 
which boosts the CM-system of the W bosons and causes a migration of events from 
regions with large cross section in the CM system to regions with small cross section in the 
laboratory system. The relative corrections are very large there, because the lowest-order 
cross section is very small. Owing to the small cross section, this region of phase space is 
rather unimportant. Similar but less dramatic boost effects are also visible in the following 
distributions. The production-angle distribution depends on the recombination scheme 
only weakly, as expected. An increase of the recombination cut leads to a redistribution 
of events from the backward to the forward direction. Since there is more phase space for 
hard photons for higher scattering energies, the effects of a full recombination are more 
pronounced at yfs = 500 GeV than at LEP2 energies, discussed in Ref. [ |I0| . This feature 
can also be observed in the following distributions. 

The distribution of events in the angle between the W + and the outgoing /z + is pre- 
sented in Figure |2|. In this case we find a noticeable dependence on the recombination 
mass M rec for large decay angles, where the cross section is, however, very small. This 
originates from the fact that the recombination of a fermion with a photon parallel to 
this fermion decreases the angle between the fermion and the W boson from which the 
fermion results. 

The distribution of events in the energy of the outgoing /i + is depicted in Figure ||. 
In the on-shell approximation it would be restricted between 6.6 GeV < < 243.4 GeV. 
Again the radiative corrections tend to smooth the distribution. The recombination of 
a photon with the muon increases the muon energy. Consequently, an increase of the 
recombination mass M rec shifts the distribution to larger muon energies, as can be seen 
in the relative corrections. 
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Figure 3: Muon-energy distribution for e + e — > z^/z+du and \fs = 500 GeV 
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Finally, in Figures |] and |5] we show the distributions of events in the invariant masses 
of the final-state lepton pair, M^, and of the final-state quark pair, Md u . The depen- 
dence on the recombination cut is sizeable everywhere. The results for the invariant-mass 
distributions can be understood as follows. For small recombination cuts, in most of the 
events the W bosons are defined from the decay fermions only. If a photon is emitted 
from the decay fermions and not recombined, the invariant mass of the fermions is smaller 
than the one of the decaying W boson. This leads to an enhancement of the distribution 
for invariant masses below the W resonance. This effect becomes smaller with increas- 
ing recombination mass. The enhancement is proportional to the squared charges of the 
final-state fermion, i.e. it is largest for the leptonic invariant mass. On the other hand, if 
the recombination mass gets large, the probability increases that the recombined fermion 
momenta receive contributions from photons that are radiated during the W-production 
subprocess or from the decay fermions of the other W boson. This leads to positive cor- 
rections above the considered W resonance. The effect is larger for the hadronic invariant 
mass since in this case, two decay fermions (the two quarks) can be combined with the 
photon. The effect of the squared charges of the final-state fermions is marginal in this 
case because the contribution of initial-state fermions dominates. 

In contrast to LEP2 energies [ ITU there is a significant difference between the loose 
recombination cut M rec = 25 GeV and the full recombination M rec = 500 GeV. As al- 
ready explained in Ref . [ [HI] , the recombination of a photon that forms an invariant mass 
of at least 25 GeV with charged fermions shifts the events by at least 3.8 GeV for the 
leptonically-decaying W boson and 7.4 GeV for the hadronically-decaying W boson. Such 
a shift leads to an increase of the corrections only outside the range shown in Figures |] 
and H However, in contrast to the LEP2 energy, the phase space for hard photons is very 
large at yjs = 500 GeV, leading to a global significant decrease of the W line shape near 
the resonance if the recombination cut M rPC increases from 25 GeV to 500 GeV. 



6 



da 



AM, 



fit/ fX 



fb 
GeV 



i — i — i — r 

Born 
5 GeV 
25 GeV 
500 GeV 




i I I I I I I I I l_ 

75 76 77 78 79 80 81 82 83 84 85 



[GeV] 



5 [% 
20 



15 



10 - 



5 - 



- 



"i — i — i — i — i — i — i — i — r 



i _ 



J I I I I I I I L 



75 76 77 78 79 80 81 82 83 84 85 
M_ [GeV] 



Figure 4: Invariant-mass distribution of the lepton pair for e + e — > z/^/i+du and y/s 
500 GeV 
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Table 2: Results for Mw,fit for a Breit-Wigner fit to the invariant-mass distributions 
shown in Figures f| and |5|, using the fit ranges 78.3 GeV < M < 82.3 GeV (upper values) 
and 76.3 GeV < M < 84.3 GeV (lower values) 



The distortion of the W invariant-mass distributions is of particular interest for the 
reconstruction of the W-boson mass M\y from the decay products. In order to illustrate 
the impact of the corrections on the determination of Mw, we fit the simple Breit-Wigner 
distribution 

/ da \ const. 

VdMV = W ~ Mlw)~ 2 + ( } 

to the W line shapes shown in Figures £| and ^ treating M Wjfit and r W fit (as well as the 
normalization constant) as free fit parameters. We determine the fitted W-boson masses 
from the predictions resulting from the lowest-order CC03 diagrams, M^ ^' 0003 , from the 
complete lowest-order diagrams, M 6 /^, and from the fully corrected predictions, M^ r g t . 
In addition we give the mass shift resulting from the corrections AM^j t = M"j t — M^et- 
The results of these fits, which are contained in Table show that the fitted W-boson 
mass changes at the order of some 10 MeV if the corrections are included. The mass 
shift crucially depends on the recombination procedure. From the discussion of the line- 
shape distortion above it is clear that this mass shift is more positive if more photons are 
recombined. Moreover, the differences in M^j t between the recombination cuts M rec = 
25 GeV and 500 GeV are only marginal, since the corresponding line-shape distortions are 
practically the same. The results also illustrate that the fit results vary at the order of 
some 10 MeV for different fit ranges in M. 

In summary, we have applied the event generator RacoonWW for e + e~ — > WW — > 
4/ (+7), which includes the complete 0(a) electroweak corrections in double-pole approx- 
imation, for the linear-collider centre-of-mass energy 500 GeV. With this generator we 
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have calculated the total cross sections and various distributions of experimental rele- 
vance. For angular and energy distributions we find corrections up to several 10%, i.e. 
they are much larger than at LEP2. Since the bulk of these corrections is due to a redis- 
tribution of events by kinematical effects induced by hard-photon radiation, corrections 
of even higher order are expected to be less drastic. For invariant-mass distributions of 
the W bosons, the size of the corrections is similar as at LEP2. The detailed numerical 
discussion demonstrates that the issue of photon recombination is of great importance in 
view of the aimed accuracy of 15 MeV in the kinematical reconstruction of Mw at TESLA 
HOI. 
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